Microscopic Spectrum of the Wilson Dirac Operator 



O 

(N 

Oh! 
CD 

o 

(N 



CD 

> 
m 

On 
(N 

o 

O 



X 



P.H. Damgaard, 1 K. Splittorff, 2 and J.J.M. Verbaarschot 3 

'Niels Bohr International Academy, Niels Bohr Institute, 
Blegdamsvej 17, DK-2100, Copenhagen 0, Denmark 
Niels Bohr Institute, Blegdamsvej 17, DK-2100, Copenhagen 0, Denmark 
3 Department of Physics and Astronomy, SUNY, Stony Brook, New York 11794, USA 

(Dated: September 21, 2010) 

We calculate the leading contribution to the spectral density of the Wilson Dirac operator using 
chiral perturbation theory where volume and lattice spacing corrections are given by universal scaling 
functions. We find analytical expressions for the spectral density on the scale of the average level 
spacing, and introduce a chiral Random Matrix Theory that reproduces these results. Our work 
opens up a novel approach to the infinite volume limit of lattice gauge theory at finite lattice spacing 
and new ways to extract coefficients of Wilson chiral perturbation theory. 



Introduction. Spectral gaps and their suppression by dis- 
order are essential for a variety of physical phenomena. 
States that intrude into the band gap, so called Lifshitz 
tail states, affect the conductivity of semiconductors [l|, 
they may lead to gapless superconductivity in supercon- 
ductors with magnetic impurities Q , and they may show 
universal fluctuations given by Random Matrix Theory 
(see 0|). Here we analyze the spectrum of the Wil- 
son Dirac operator of lattice Quantum Chromodynamics 
(QCD). In the continuum limit, the Hermitian Wilson 
Dirac operator has a gap equal to twice the quark mass. 
At finite lattice spacing, eigenvalues of tail states intrude 
into the gap. When eigenvalues approach the center of 
the gap, it becomes increasingly difficult to invert the 
Wilson Dirac operator. As a consequence, such tail states 
can potentially obstruct lattice simulations. It is there- 
fore of importance to have an analytical understanding 
of the properties of these states. 



Our results rely on two approaches, chiral Random 
Matrix Theory and chiral Lagrangians for the pseudo- 
Nambu-Goldstone sector of QCD. The relation between 
chiral Random Matrix Theory and the Dirac operator 
in theories with spontaneously broken chiral symmetries 
Q has led to a new understanding of the chiral limit 
of strongly coupled gauge theories. The Random Ma- 
trix Theory results are universal [f| and are equivalent 
Q to what is obtained from a chiral Lagrangian in the 
microscopic domain or e-regime Q- This gives a finite- 
volume scaling theory for spectral correlation functions 
as well as individual eigenvalue distributions of the con- 
tinuum Dirac operator at fixed topological charge v. In 
lattice QCD it has become standard to utilize these re- 
sults to obtain physical observables from simulations at 
finite four- volume V. There has for long been a desire 
to obtain analogous results for Wilson fermions at finite 
lattice spacing a. Here we present a solution to this prob- 
lem. 



We denote the Wilson Dirac operator by D = Dw - 



Below we make use only of its block structure given by 

with = A and B< = B, whereas W does not have 
additional symmetry properties. The Wilson Dirac oper- 
ator is anti-Hcrmitian in the continuum limit a — > and 
the corresponding eigenvalues of Dw are complex away 
from the continuum limit. 

Although non-Hermitian, the Wilson Dirac operator 
satisfies 75-Hcrmiticity 



75-D75 



(2) 



Instead of the Wilson Dirac operator itself, it is therefore 
often more convenient to work with the Hermitian Dirac 
operator D5 = J5D. At zero lattice spacing, the spec- 
trum of -D5 has a gap around the origin of width 2m. At 
non-zero lattice spacing a, states intrude inside the gap 
and for sufficiently large lattice spacing the gap closes. 
Then one enters what is known as the Aoki phase [f|. 
It is reminiscent of the Gorkov Hamiltonian for super- 
conductors, where magnetic impurities (see @, H|) pl a y 
the role of the diagonal blocks in Eq. (p} and the Aoki 
phase corresponds to gapless superconductivity. A first 
order scenario where the condensate jumps as a function 
of m has been suggested Q and support for this has been 
found on the lattice ■ 

The spectral density p^(x) of D5, evaluated at x = 0, is 
an order parameter ll| for the onset of the Aoki phase. 
Discretization effects in the spectrum of the Wilson Dirac 
operator were analyzed by means of chiral perturbation 
theory in [l2| • What is new here is that we obtain an ex- 
act analytical description in the microscopic scaling limit 
and show in detail the transition to the Aoki phase. This 
opens a novel analytical approach to the infinite-volume 
limit of lattice gauge theory at finite lattice spacing and 
offers new ways to measure the leading coefficients of Wil- 
son chiral perturbation theory. Understanding the distri- 
butions of the low-lying eigenvalues of the Wilson Dirac 
operator is also crucial for establishing a stable domain 
for numerical simulations [l3| . 
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Chiral Lagrangian. The leading-order terms of the chiral 
Lagrangian for Wilson fermions have been listed in Q. 
It is a double expansion: the continuum ordering for chi- 
ral perturbation theory and an expansion in the lattice 
spacing a. The corresponding chiral Lagrangian coin- 
cides with the continuum Lagrangian with shifted mass 
plus terms starting at order a 2 . It is convenient to in- 
troduce a source for ipjsip, which we denote by z. Here, 
we shall focus on the microscopic domain where mV, zV 
and a 2 V are kept fixed in the infinite- volume limit. Dif- 



ferent counting rules are possible |14| , but the present one 
is most useful for elucidating the effects of finite lattice 
spacing on the low-lying Dirac eigenvalues. The leading 
contribution to the finite- volume QCD partition function 
then reduces to a unitary matrix integral, which, up to 
a few constants, is determined by symmetry arguments. 
We decompose this partition function as Z^ f = Y],, Z^^ 
with 

Z v N {m,z;a) = [ dUdet v Ue s M (3) 

Ju(N f ) 

where the action S[U] for degenerate quark masses is 

S = y£FTr([/ + [/ t ) + ^£VTr([/-C/ t ) (4) 

-a 2 VW 6 [Tr (U + U^)} 2 - a 2 VW 7 [Tr (U - U^)] 2 

~a 2 VW s Ti(U 2 + U^). 

Below we will demonstrate that in the microscopic do- 
main Zpj corresponds to ensembles of gauge field config- 
urations with v real modes of D\y . The a 2 -terms are de- 
termined by invariance arguments 0], and E, We, W7 and 
Ws are the low-energy constants of 0(a 2 ) Wilson ChPT. 
The two terms corresponding to Wq and W7 are expected 
to be suppressed in the large- N c limit [TBI, and we shall 
for simplicity i gno re them here. The potential impact 
of these terms [l6| can be studied at the expense of a 
slightly more cumbersome analysis. The leading finite- 
volume partition function Zf^ then only depends on the 
microscopic scaling variables to = mSV, z = zHV, and 
a = a\/WsV which will be kept fixed for V — > 00. The 
sign of Ws will be discussed below. 

The Generating Function. A generating function for 
spectral correlation functions is given by an average of 
ratios of determinants. Because of the inverse determi- 
nants, it has an extended graded flavor symmetry. The 
graded generating function for spectral correlations of D§ 
is 



Z v w {M,Z;a) = JdU Sdet(U) v 



(5) 



x i^Str(MlU-U' 1 ])+i^Str(Z[U+U' 1 ])+aStr(U 2 + U' 2 ) 

where M. = diag(mi . . . rhk+i) and Z = diag(ii . . . Zk+l)- 
This graded partition function differs in a subtle way 
from the one introduced in 11211 . As discussed in |6f, the 
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FIG. 1: The microscopic spectrum of D5 for rh = 3, v = 
and a = 0, 0.03, and 0.250. The v — spectrum is reflection 
symmetric about x — 0. 



integration manifold for non-perturbative computations 
is non-compact for the bosonic sector. While the action 
(HJ) and the action introduced in [l2| break the flavor 
symmetries in exactly the same way, only the former one 
is consistent with the convergence requirements of the 
graded integral for Ws > 0. For pcrturbativc calcula- 
tions the convergence requirements are immaterial. Here 
and below we focus mainly on the quenched case, corre- 
sponding to integration manifold GI(\\1)/U{1). The gen- 
eralization to an arbitrary number of flavors is straight- 
forward, and we expect that the underlying integrability 
structure will lead to a full analytical solution just as in 
the a = case 17 1 . 

The Microscopic Spectrum of D$. For a = 0, the micro- 
scopic spectral density of D5 follows from the expression 
for the microscopic spectral density of D through 



p${x > rh, m; a = 0) = 



\J x 2 — rh 2 



p u {yjx 2 -rh 2 ) . (6) 



To obtain the spectral density of D 5 for o ^ 0, we eval- 
uate the resolvent 



G u (z,rh]a) = lim — Z^ijh, rh, z, z ; a) (7) 

2'— >i dz. 



and find 



G v {z, rh; a) 



00 ds ^ L cos{ 0) e Sf+S be (;W-s), 

-00 J-ir 27r 2 



x (—to sin(0) + irh sinh(s) + ii cos(0) + iz cosh(s) (8) 
+4a 2 [cos(26») + cosh(2s) + (e l0+s + e" 49 " 8 )] + l) . 

Here Sf = —rhs'm(8) + izcos(8) + 2d 2 cos(26>) and S& = 
—irh sinh(s) — iz cosh(s) — 2a 2 cosh(2s) . As is well known, 
the resolvent is defined only up to ultraviolet subtrac- 
tions in the underlying theory. The microscopic quenched 
spectral density, 



/9c (x, rh; a) = — Im[G l/ (x, to; &))], (9) 

7T 
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divergent. However, the graded partition function 



FIG. 2: The microscopic spectrum of D5 for rh = 3, v 
and a = 0.125, a — 0.250 and a = 0.500 respectively. 



is, however, uncontaminated by these ultraviolet pieces. 
Plots of are shown in Fig.[T]for v = 0, and in Fig.[2]for 
v = 2. The eigenvalues that converge towards the end- 
points of the spectrum at sign(z/)m in the a — >• limit are 
clearly visible. The sum over v of the spectral density, 
for which the continuum limit has been established rig- 
orously [lH , can be evaluated in a straightforward way. 
Random Matrix Theory. An efficient alternative way to 
extend the above results to all spectral correlation func- 
tions and individual eigenvalue distributions is to con- 
struct a chiral Random Matrix Theory that is equiva- 
lent to the chiral Lagrangian in the same scaling regime. 
This the case if the chiral Random Matrix Theory has 
the same global symmetries and transformation proper- 
ties as the QCD partition function with the Wilson Dirac 
operator. The chiral Random Matrix Theory is 

Z V N = I dAdBdW det Nf {D w +m + z 75 ) P(D W ), (10) 



where Dw is of the same block form as ((T|) and the inte- 
gration is over the real and imaginary parts of the matrix 
elements of the Hermitian nxn matrix, A, the Hermitian 
(n + v) x (n + v) matrix B and the complex n x (n + v) 
matrix W. This Dw has \v\ real eigenvalues. We have 
added tildes to stress that this is a zero-dimensional ma- 
trix integral with parameters rh and z instead of m and 
z. In the universal scaling limit there is a one-to-one cor- 
respondence between the two pairs, just as in the a = 
case. The precise form of the distribution of the matrix 
elements P(Dw) is not important on account of univer- 
sality. The partition function §5§ (with W$ = Wj = 0) is 
recovered in the microscopic scaling limit. For a Gaus- 
sian distribution, this can be shown by a simple explicit 
calculation. It also follows that Wg > 0. This is a con- 
sequence of the 75-Hermiticity: Changing the sign of Wg 
is equivalent to a — > ia, violating 75-Hermiticity of Dw- 
This suggests that the Hermiticity properties of the Dirac 
operator can restrict the coefficients of the effective La- 
grangian. In fact, with Wg < the integrals in ([5]) are 



Z^MiZ-a) = JdUSdet(uy (11) 

,^Str(M[U+U' 1 ]) + ^StiiZ[U-U' 1 ])~a 2 Str{U 2 + U^ 2 ) 



xe 2 " 



is now convergent. Repeating the steps leading to 
with this convergent integral, we find a resolvent for an 
operator that, unlike D5, is not Hermitian. We believe 



that the absence of solutions in the p-regime [12| for Wg < 
has the same origin. 

The ensemble with the structure of the matrix Dw 
belongs to one of the classes in the non-Hermitian clas- 
sification of [l9[ (the 75-Hermiticity is there referred to 
as Q" s y mme try) . In the microscopic scaling limit, the 
partition function (|10j) has the determinantal structure 



Z v Nj (fh,z;a) = det[Z^ 1 J (?7i, z; a)] iJ=1 ... Nf (12) 



where 



I e i6»i/ e ?ricos(e)+izsin(e)-2o 2 cos(20) 



This form suggests that the partition function is a r- 
function of an integrablc system of Toda type and that 
an eigenvalue representation can be obtained. 

The simplest quantity to compute from (|12j) is the chi- 
ral condensate, S(m) = dm log Z. For v = there is 
a striking similarity to the chiral condensate for QCD 
at non-zero isospin chemical potential /x. The conden- 
sate is constant for large rh and drops roughly linearly 
to zero inside a well defined region (the Aoki phase and 
the pion condensed phase, respectively). This is not ac- 
cidental: The microscopic spectrum of Dw at a =/= 
forms a thin strip along the imaginary axis just as the 
continuum Dirac operator does at n 7^ 0. In both cases, 
the chiral condensate can be interpreted as the electric 
field at m of point charges at the position of the eigen- 
values. Also the convergence requirements of the graded 
partition function d5j have direct analogues at nonzero 
chemical potential [20(. 

The Density of Real Modes. The analytical result for 
the quenched average spectral density of the real eigen- 
values, denoted by p^ rcal , also follows from the generat- 
ing function §5§ (a derivation will be given in [21(). A 
plot of p^ rcal for v = 4 versus ( = (Y^V is shown in 
Fig- El The real eigenvalues, of Dw repel each other 
and the \v\ = 4 real modes are clearly visible. We have 
checked that the distribution is in agreement with the 
chiral Random Matrix Theory (fTOj) . This illustrates that 
in the label v introduced in ([3]) corresponds to the num- 
ber of real eigenvalues. For large v we find that Pwroai 
approaches a semi-circle. We suggest that analyzing just 
these real eigenmodes of Dw may provide a new useful 
tool in lattice QCD. 
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FIG. 3: The quenched density of the real eigenvalues of Dw- 



tial to this study is an analysis of the gap of the Wilson 
Dirac operator at finite mass. Lattice QCD simulations 
depend crucially on control of this gap and its variation 
as a function of the lattice spacing. As we have stressed, 
our results are not only important for understanding the 
finite lattice spacing effects of Wilson fermions near the 
chiral limit, but may have interesting applications to tail 
states in condensed matter systems. 
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Distribution of Tail States. For \x — rh\/a 2 3> 1 and 
8a 2 <C 1 the tail of the spectral density inside the gap 
follows from a saddle point analysis. For x > we find 

P5 {x) ~cxp[-(x-m) 2 /16a 2 ], (14) 

and a similar result for x < 0. This result applies to the 
tail of the blue and (marginally) of the red curve in Fig.[T] 
Reinstating physical parameters we find that the width 
parameter, a, in (|14[) is given by a 2 = 8a 2 Ws/(VT, 2 ) so 
that, in the microscopic domain and sufficiently small a, 
a scales with \l\fV . Such a scaling has been observed 
for N f = 2 in []j,[22|. 

Tail states may also be studied in their own right and 
for applications in condensed matter physics. In particu- 
lar, in the thermodynamic limit, rh,x,a 2 3> 1, the aver- 
age level density of D$ can be obtained by a saddle point 
analysis. For 8a 2 /in < 1 it vanishes inside [— x c ,x c ] with 
x c given by x c = 8a 2 [(m/8a 2 ) 2 ^ 3 — l] 3 / 2 . It has exactly 
the same form as for superconductors with magnetic im- 
purities 0. In the scaling limit where V 2 ^ 3 (x c — x) is 
kept fixed, the spectral density can be computed inside 
the gap by a saddle point approximation of ([5]), and it 
agrees with universal Random Matrix Theory results for 
the so-called soft edge. Such universal behaviour has also 
been found in condensed matter systems @, Q . 
Conclusions. Using a graded chiral Lagrangian for Wil- 
son fermions at finite lattice spacings, we have obtained 
an analytical form for the Wilson Dirac spectrum at fixed 
number, v, of real eigenvalues. These results, and their 
extensions to dynamical fermions, should be useful for 
lattice simulations at finite volume. We have shown how 
the leading low-energy constant for Wilson fermions, W$, 
can be extracted from lattice spectra of the Wilson Dirac 
operator in the e-regime. 

The problem can also be reformulated in terms of a 
new chiral Random Matrix Theory that describes spec- 
tral correlation functions of the Wilson Dirac operator 
in the appropriate scaling regime. These results open 
a new domain of Random Matrix Theory where chiral 
ensembles merge with Wigner-Dyson ensembles. Esscn- 



[1] I.M. Lifshitz, Sov. Phys. Usp. 7, 549 (1965). 

[2] A. Lamacraft and B. D. Simons, Phys. Rev. Lett. 85, 
4783 (2000); Phys. Rev. B 64, 014514 (2001). 

[3] C.W.J. Beenakker, Lect. Notes in Phys. 667 131 (2005). 

[4] E.V. Shuryak and J.J.M. Verbaarschot, Nucl. Phys. 
A560, 306 (1993); J.J.M. Verbaarschot, Phys. Rev. Lett. 
72, 2531 (1994). 

[5] G. Akemann, P. H. Damgaard, U. Magnea and S. Nishi- 
gaki, Nucl. Phys. B 487, 721 (1997). 

[6] P. H. Damgaard, J. C. Osborn, D. Toublan and 
J. J. M. Verbaarschot, Nucl. Phys. B 547, 305 (1999). 

[7] J. Gasser and H. Leutwyler, Phys. Lett. B 188, 477 
(1987); H. Leutwyler and A. V. Smilga, Phys. Rev. D 
46, 5607 (1992). 

[8] S. Aoki, Phys. Rev. D 30, 2653 (1984). 

[9] S. R. Sharpe and R. L. Singleton, Phys. Rev. D 58, 
074501 (1998); Rupak and N. Shoresh, Phys. Rev. D 
66, 054503 (2002); O. Bar, G. Rupak and N. Shoresh, 
Phys. Rev. D 70, 034508 (2004); S. R. Sharpe and 
J. M. S. Wu, Phys. Rev. D 70, 094029 (2004); M. Goiter- 
man, S. R. Sharpe and R. L. Singleton, Phys. Rev. D 71, 
094503 (2005). 

[10] F. Farchioni et at, Eur. Phys. J. C 39, 421 (2005); 

C. Michael and C. Urbach, PoS LAT2007, 122 (2007); 

P. Boucaud et al., Comput. Phys. Com. 179, 695 (2008). 
[11] K. M. Bitar, U. M. Heller and R. Narayanan, Phys. Lett. 

B 418, 167 (1998). 
[12] S. R. Sharpe, Phys. Rev. D 74, 014512 (2006). 
[13] L. Del Debbio, L. Giusti, M. Luscher, R. Petronzio and 

N. Tantalo, JHEP 0602, 011 (2006). 
[14] A. Shindler, Phys. Lett. B 672, 82 (2009) O. Bar, 

S. Necco and S. Schaefer, JHEP 0903, 006 (2009). 
[15] R. Kaiser, H. Leutwyler, Eur. Phys. J. C 17, 623 (2000). 
[16] S. R. Sharpe, Phys. Rev. D 79, 054503 (2009). 
[17] K. Splittorff and J. J. M. Verbaarschot, Phys. Rev. Lett. 

90, 041601 (2003). 
[18] L. Giusti a nd M. Luscher, JH EP 0903, 013 (2009). 
[19] U. Magnea. larXiv:0707.0418V 2 [math-ph]. 
[20] K. Splittorff and J. J. M. Verbaarschot, Nucl. Phys. B 

683, 467 (2004); Nucl. Phys. B 757, 259 (2006). 
[21] G. Akemann, P.H. Damgaard, K. Splittorff and J.J.M. 

Verbaarschot, to appear. 
[22] L. Del Debbio, L. Giusti, M. Luscher, R. Petronzio and 

N. Tantalo, JHEP 0702, 082 (2007). 



